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Abstract. The formation of correlations due to collisions in an interacting nucleonic system is investigated. 
Results from one-time kinetic equations are compared with the Kadanoff and Baym two-time equation with 
collisions included in Born approximation. A reasonable agreement is found for a proposed approximation 
of the memory effects by a finite duration of collisions. This form of collision integral is in agreement 
with intuitive estimates from Fermi's golden rule. The formation of correlations and the build up time is 
calculated analytically for the high temperature and the low temperature limit. Different approximate ex- 
pressions are compared with the numerical results. We present analytically the time dependent interaction 
energy and the formation time for Gaufi- and Yukawa type of potentials. 
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1 Introduction 

The Boltzmann transport equation has played a very im- 
portant role in the development of non-equilibrium statis- 
tical mechanics. This microscopic equation describes the 
time-evolution of a distribution-function in phase-space 
and has also provided a connection with macroscopic hy- 
drodynamic equations by a moment expansion of the mo- 
mentum. Important applications are for example the well- 
known Chapman-Enskog calculations of transport coeffi- 
cients. In later developments the Markovian Boltzmann- 
equation has been extended to include memory and correlati 
effects in the collision-integral and there are a large num- 
ber of publications concerning such improvements. These 
classical kinetic equations describe the time-evolution of 
a one-time distribution function /(r, p, t). 

Meanwhile, a quantum two-time theory for the time- 
evolution of real time Green's functions G(r, p, t, t') has 
been developed using the Schwinger-Keldysh formalism. 
The quantum image of the classical Boltzmann equation 
is usually referred to as the Kadanoff-Baym (KB) equa- 
tions J|] . These equations have often been considered too 
complicated to solve numerically in the past. However, 
several numerical applications exist now. The Kadanoff- 
Baym equations have also played an important role in 
the improvements of the Boltzmann equation especially 
by using the Generalised Kadanoff Baym Ansatz (GKB) 
of Lipavsky et al . This ansatz allows a reduction of the 
two-time formalism to a formally simpler one-time formal- 
ism e.g. the Boltzmann equation. The time off-diagonal 
Green's function elements are related by GKB to the time- 
diagonal by the spectral-functions. By various approxima- 



tions of the spectral-functions various one-time approxi- 
mations of the two-time equations can be obtained. (See 
e.g. |) 

These kinetic equations describe different relaxation 
stages. During the very fast first stage, correlations im- 
posed by the initial preparation of the system are decay- 
ing J|,||. These are contained in off-shell or dephasing 
processes described by two-time propagators. During this 
stage of relaxation the quasiparticle picture is established 
[pin]] . After this very fast process the second state develops 
cfuring which the one-particle distribution relaxes towards 
the equilibrium value with a relaxation time t ic \. First 
the momentum anisotropy relaxes by small angle scatter- 
ing events and then the energetic degrees of freedom relax. 
During this relaxation state the virial corrections are es- 
tablished and can be consistently described by a nonlocal 
Boltzmann kinetic equation |||l0|. The time of the first 
stage t c is mostly shorter than the relaxation time T re i of 
one particle distributions which is entirely determined by 
the collision process. We will focus on the first stage which 
is related to the formation of correlations. 



The formation of correlations is connected with an in- 
crease of the kinetic energy or equivalently the build up of 
correlation energy. This is due to rearrangement processes 
which let decay higher order correlation functions until 
only the one - particle distribution function relaxes. Be- 
cause the correlation energy is a two - particle observable 
we expect that the relaxation of higher order correlations 
can be observed best within this quantity. Of course, the 
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total energy of the system is conserved 



(i) 



Chapter IV summarizes our results and we discuss 
some aspects regarding the correlation time. The appen- 
dices show some important relations necessary for our an- 
alytic calculations. 



which means that the kinetic energy increases on cost of 
the correlation energy E COTI (t). We will observe a transfor- 
mation of correlation into kinetic energy. This process sat- 
urates on the end of the first stage of relaxation. It is more 
convenient to calculate the kinetic energy than the corre- 
lation energy because the kinetic energy is a one-particle 
observable. Consequently, the time dependence of the ki- 
netic energy will be investigated within the kinetic theory. 
This can only be accomplished if we employ a kinetic equa- 
tion which leads to the total energy conservation (|l|) . It is 
immediately obvious that the ordinary Boltzmann equa- 
tion cannot be appropriate for this purpose because the 
kinetic energy is in this case an invariant of the collision 
integral and constant in time. Imposing the conservation 
of the form (|l|) we have to consider non-Markovian ki- 
netic equations JlT| , which account for the formation of 
two particle correlations. 

Within these kinetic equations the collision integral 
is an expression of the two-particle correlations. While 
the one-particle distribution remains almost unchanged 
during the first stage of relaxation, the two-particle cor- 
relations relaxes. Consequently the one- particle spectral 
function is changing. The latter one is responsible for the 
dephasing and therefore formation of correlations. Ana- 
lytical expressions for the time dependence of the kinetic 
and correlation energy are obtained in this paper by con- 
sidering explicitly this dephasing process. 

We start from a kinetic equation appropriate for short 
time scale in Born approximation. It contains the full 
memory-effect but no damping i.e. no explicit width of 
the spectral function, because quasiparticles are not yet 
formed on this time scales. In Chapter II we give an overview 
of the gradient approximation with emphasis on energy- 
conservation and correlation energy. In appendix [b| we 
discuss the limit of complete collisions and the weaken- 
ing of initial correlations. In appendix ^ we calculate the 
equilibrium value of the correlation energy for high and 
low temperature limits analytically using Gaussian and 
Yukawa type interactions. 

On the very first time scale we can neglect retardation 
effects in the one particle distribution function, but we 
have to keep into account off-shell properties of the colli- 
sion integral. Therefore we use the finite duration approx- 
imation in Chapter III. It leads to the correct equilibrium 
value and is intuitively clear from Fermi's Golden rule. It 
is compared numerically with the KB results. These cal- 
culations also bring to the attention the correlation-time 
i.e. the time for build-up of correlations. 

From the observation that the time-variation of the 
distribution-functions can be neglected in the first stage 
of relaxation we obtain an analytic expression for the time 
dependent formation of correlations. Especially we give 
analytical results for the formation time of correlations in 
a high and in a low temperature limit. Comparisons are 
made with numerical calculations. 



2 Correlation Energy in gradient expansion 

The kinetic equation in Born approximation for spatial 
homogeneous media including complete time convolution 
(memory effect) but no damping is called Levinson equa- 
tion and reads |fiM3M[L5| 
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with / = 1 — /, the free particle dispersion E = p 2 /2m 
and the spin-isospin degeneracy Si,s 2 . The distribution 
functions are normalized to the density as s J p^fp f(p) = 
n. For the sake of simplicity we have omitted the Hartree 
and Fock contribution. Since we discuss the correlation 
energy the meanficld contribution is just additive. The 
dispersion E{p) in the collision integral is modified by the 
Fock term, but we use this effect only in an approximative 
way by understanding m as effective mass. 

The Boltzmann collision integral is obtained from equa- 
tion (H) if: (i) One neglects the time retardation in the dis- 
tribution functions, i.e. the memory effects and (ii) The 
finite initial time to is set equal to — oo corresponding to 
what is usually referred to as the limit of complete colli- 
sions. The memory effect is condensed in the explicit re- 
tardation of the distribution function. This would lead to 
gradient contributions to the kinetic equation which can 
be shown to be responsible for the formation of high ener- 
getic tails in the distribution function [jl6,17 . This effect 
will be established on the second stage of relaxation. 

The second effect is contained in the energy broaden- 
ing or off-shell behavior in (0). This is exclusively related 
to the spectral properties of the one-particle propagator 
and therefore determined by the relaxation of two-particle 
correlation. Since we are studying the very short time re- 
gion after the initial disturbance we can separate the one- 
particle and two-particle relaxation. On this time scale the 
memory in the distribution functions can be neglected but 
we will keep the spectral relaxation implicit in the off-shell 
cos-function of (||) . This effect is the most relevant one for 
obtaining the time evolution of the interaction (or corre- 
lation) energy and therefore energy conservation. 

In the following discussion we shall only be concerned 
with the time integration. Therefore we introduce the short 
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hand notation of equation 



Ft f(pi)= M to drcos 



of the collision integral directly from Eq. 
to 



This leads 



AE(t - t) 



F(t), (3) 



di 



E co „ = J2<I n (t)F^(t) > 



where 



F(t) =2V(\ vx~ P'x \) 2 5{PI+P2-P 1 -P 2 ) 

x(f(p'l,T)f( P ' 2 ,T)f( Pl ,T)f(p 2 ,T) 
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dt't' n AE cos 



AEt 
h ' 



(9) 



(4) 



and the 9-dimensional momentum-integration is suppressed 
in Eq. (||) and in the following. 

From this equation one derives balance equations by 
integration over momentum p\. The first two moments, i.e. 
the density and total linear momentum, are conserved. For 
the Markovian Boltzmann equation the kinetic energy is 
conserved, while potential energy is zero. In the present 
case including the memory effect one finds ftl| 



!«ji> +E „„) = o, 

where the correlation energy E cm T- is given by 



(5) 



Note the difference between the two expansions @ and 
(||). For example, the zero order term in Eq.(||) does con- 
tain not only the zero order term but also part of the first 
order term of Eq.(^). This is understandable, because the 
collision integral determines the time derivative (|^) of the 
correlation energy. One has to expand the collision integral 
one step further in order to obtain the correlation energy 
(0) up to a specific level of gradient expansion. This is 
a quite general observation for any order of gradient ap- 
proximation. 

Comparing the two gradient expansions (^) and (|^) we 
establish a relation between I n (t) and V n (t) 



1 / ,t "*° f AEt 

E co „(t) ~ E COTT (to) = -|^( / dT sin ( — — ) F ( f - T ))- 
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(10) 



(6) or inversely 



Here <> indicates the integration over p\. We like to 
point out that we have neglected initial correlations in 
the kinetic equation (^) in agreement with the studied 
sudden switching approximation. Consequently, E coir (to) 
describes only possible constant background correlations 
not formed by the binary collisions. Expanding F(t — r) 
around t one obtains a gradient expansion series for the 
interaction energy that reads 

oo 

E COII (t) ~ ^corr(O) = < V n (t)F {n \t) > 



V n (t) = J dt'(I n (t')-V n ^(t')) 

to 

i -. , AE(t-t ) 

v (t) = J dfw) = i cos 2 E 

to 



(11) 



n=0 



v n (t) = - 



1 (-1)" 



dt't' n sin 



AEt' 



The long time limit of the different gradient approxima- 
tions of the kinetic equation are presented in appendix [b| 
and is found to be unique. The limit of complete colli- 
sions to — > — oo and the connected problem of weakening 
of initial correlations are discussed there. 



( 7 ) 3 The formation of correlations 



where the n-th time derivative of F(t) is given by (t) 
Taking the time-derivative of (m) one finds 



dl Eco " = ^ 

n=0 



< V£(i)F< n >(i) + V n (t)F {n+1 \t) > 



(8) 



On the other hand one can express the time derivative 
of the interaction energy in terms of a gradient expansion 



To lowest order the gradient expansion (Q) , the correlation 
energy is 



I CO g AE{t-tg) _ ^ 

E COII (t) - Scorr(O) = -( 



F<®(t)). 



(12) 



Retaining only the first term for this correlation energy 
is equivalent to an approximation of the non-Markovian 
collision integral (Q) where we neglect the time dependence 
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of the distribution functions while keeping the finite initial 
time to- This approximation gives instead of (|^) 

„ . AE(t-t ) 



2 f dp 2 dp[dp 2 



V(\pi-p[ \) 2 6(pi+ P 2-p' 1 -p' 2 ) 



k J (27T?i) 6 

sin (Ei + E 2 -E[- E' 2 ){t - t )/% 
Ei + E% — E[ — E' 2 

X [f(p'i,t)f(jp' 2 ,t)f(#x,t)f(jp2,t) 
- f(pl,t)f(p2,t)f(p[,t)f(p 2 ,t)]. 



(13) 



Using the same steps which was used to derive the corre- 
lation energy E co „(t) in @ from the collision integral (|) 
one easily finds that the collision integral (13) gives the 
lowest order term of the time derivative of this correlation 
energy, i.e. (fl2]). The off-shell function in this collision in- 
tegral contains a memory of the initial state at to. This 
induces a memory in the kinetic equation in spite of the 
fact that the collision integral is formally Markovian. From 
(Q) one recognizes that the equilibrium or long time limit 
of the correlation energy is exactly give by the first order 
gradient approximation ( |l2| ) since all next orders include 
time derivatives and vanishes therefore on the long time 
scale. 

We summarize three important features of this approx- 
imate collision integral: (i) It reproduces the zero order 
gradient term of the time dependent correlation energy 
(Q). (ii) It leads to the complete expression for the equi- 
librium correlation energy (32) and (iii) it is a direct conse- 
quence of Fermi's Golden Rule, where we keep time depen- 
dent perturbation theory. Previously, this approximation 
has been considered also for electron plasmas P,|^,p"8|,|T9|]. 

It was shown in appendix |b] that the two operations 
of taking the time-derivative and the long-time limit of 
to do not commute. This is nicely illustrated by Eq. ( [l3| ) 
which gives the correct correlation energy ( |32"| ) if the limit 
to — > — oo is performed afterwards. If however the limit of 
complete collision to — > — oo (or t — > oo) is performed first 
on the kinetic equation (13) it reduces to the Markovian 
Boltzmann equation without correlation energy. 

We refer to Eq. ([l3| ) as the finite duration time approx- 
imation. It carries the most important features of the build 
up of correlations after the interactions are switched on 
in the initially uncorrelated system. This will be demon- 
strated by some numerical examples below. We shall first 
present some analytical results for high and low tempera- 
ture limit and compare them with the numerical solution 
later. To this end we assume a system consisting of two dif- 
ferent types of particles a, b with different masses m a , mt. 
As a first example we shall consider a Yukawa-type poten- 
tial of the form 



V Y {r) 



9ab c ~K,r 

r 



(14) 



where in nuclear physics applications g a b is the coupling 
and k the effective range of potential given by the inverse 



mass of interchanging mesons. In plasma physics appli- 
cations the potential ( fl4| ) represents the Debye potential 
with g a t = e a eb/eo given by the two charges and the in- 
verse screening length 



E 



Airej dn a 
£o Ma 



(15) 



where n a is the density and /z a the chemical potential of 
specie a. 

As a second example we shall use a Gaufi-type poten- 
tial 



V G (r) = V e-( r ^ 2 



(16) 



which has been used in nuclear physics applications p0| , [2l] ] 
with r) = 0.57nxT 3 and V = -453 MeV. 

We will proceed and derive analytical expressions which 
are compared with the numerical solution of ( |l2| ) . The nu- 
merical values are compared in table [l] with the solution of 
the KB equations. It is found an overall good agreement. 



3.1 Time dependent correlation energy 

We calculate the build up time of correlation, r c by in- 
specting the time derivative of the interaction energy. We 
define t c as the time at which this derivative becomes suf- 
ficiently small. This corresponds to using (||) instead of 
(0), but only with the first term according to the finite 
duration approximation of the last chapter. We have from 

& 



d_ 



E r , 



1 , . AE(t - t ) 

sin E(t)}. 



4ft 



(17) 



3.1.1 High temperature limit 

In the limit of high temperature we neglect the degener- 
acy and the equilibrium distribution takes the Maxwell- 
Boltzmann form 



f(p) = -y e -^; A 

s ml 



(18) 



We get for the Gaufi potential and b 2 = ft 2 /(2/zlV) 
the result 







—e g =y 

dt co " ^ 
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V%7Tb 4 \M 



o (^p 2 + m 2 T 2 /n 2 -i3) 

~d/3 1 8(/3 2 + l6t 2 T 2 /h 2 ) 



1/2 



l3=2/b 2 +4t 2 T 2 /h 2 



E 



Zn a n b ^' 2 ^V 2 {2^f' 2 f A^\ 2 1_ 
16T 5 / 2 \M I t 4 



+ o(t- 5 ). 



(19) 
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We see a monotonic decrease of the time derivative or 
cquivalcntly a monotonic increase of the correlation en- 
ergy to its equilibrium value. It is remarkable that the long 
time limit is entirely determined by the classical value. 
Obviously no quantum effects enter the formation of cor- 
relation in the high temperature limit. 

We can also calculate the time dependent formation 
for Yukawa like potentials. The time derivative of the cor- 
relation energy leads to [appendix (|65 ~ 



dt 



2h 



Illl 



(l + 2z 2 )e z2 (l-erf(z))-^ 

\/7T 



(20) 



where we used z = ui p \Jt 2 — it^ and the collective (plasma 

frequency to 2 — n 2 T/m, compare with (|l5|). This is the an- 
alytical quantum result of the time derivative of the for- 
mation of correlation. For the classical limit we are able to 



integrate expression (20) with respect to times and arrive 
at 



2uj p t 



(1 + 2{u p tf)e { ^ 2 (1 - erf(a; p t)) 



(.21) 



It shows that the formation of correlations is basically 
given by the inverse collective frequency oj p . Therefore 
during the first stage of relaxation the fluctuating and 
collective effects are more important than collisions. 



3.1.2 Low temperature limit 

The low temperature value is of special interest, because it 
leads to a natural definition of the build up of correlations. 



Using the same steps as in appendix (A. 2) one obtains 
from (FT 



d 



l 



v 2 



at 2 cos 



-)If 



(22) 



with abbreviation as in (Fd]) . It is now easy to perform the 
time integral in (22) to obtain the time dependence of the 
correlation energy 



E, 





with x — 2 e/ ^~ ec £ and the equilibrium correlation energy 
E 1 ™^ from ([42]) respectively. This shows that the corre- 
lation energy is built up and oscillates around the equi- 
librium value. This oscillations are damped with in 
time. 

We would like to point out here that the result for 
the formation of correlations at low temperatures (E23) is 




total 



X 



correlation 



16 



20 



time [ fm/c ] 



Fig. 1. The time dependent kinetic and correlation energy 
vs. time for a counter-flowing streams of nuclear matter from 
(|l3|). The temperatures and densities of the colliding beams 
are Ti = 10 MeV, m = n o /60 and T 2 = 5 MeV n 2 = n o /10 
and the relative momentum 1.5H/ fm, which corresponds a col- 
liding energy of 45 MeV/nucleon. The beams are starting to 
interact at time point to- One sees the build up of correlation 
energy during a correlation-time t c = 3 fm/c. The parameter 
are chosen in such a way that the system is degenerate, which is 
described by n\ 3 = 0.416 < 1. The total initial kinetic energy 
corresponds (neglecting correlation energy) to an equilibrated 
system with a temperature of T = 32 MeV. 



independent of the interaction used. Because higher or- 
der interaction described by higher order diagrams can be 
cast into a Boltzmann- like collision integral with more in- 
volved transition matrix elements we have always the same 
time dependence of (23) for binary interactions however 
with different equilibrium correlation energy E corT . 



3.2 Numerical results 



Figure [l] shows the time development of the kinetic energy 
as well as the correlation energy. The equation (|l3|) has 
been solved numerically for two initially counter-flowing 
streams of nuclear matter, where the nucleons interact via 
the Gaufi type of potential (|l6|). The initial temperatures 
and densities of the colliding beams are T\ = 10 MeV 
n\ — n o /60 and T 2 = 5 MeV 712 = n o /10 respectively. 
The relative momentum is 1.5h/ fm, which corresponds 
to a colliding energy of 45 MeV/n. We observe a build 
up of correlations during the initial 3 — 4/m/c. Total en- 
ergy is conserved and the kinetic energy is increased by 
the same amount as the correlation energy is decreased. 
Similar results has been found for a different system in 
||. This is because the system is initially prepared to be 
uncorrelated at to = 0. If time to i.e. the time when the 
system is uncorrelated is shifted to to — —00 i.e. infinite 
past, we would not observe any build up of correlations. 
The equation (|l^) would then in fact reduce to the Boltz- 
mann equation. 



,tions and energy-conservation at short time scales 




Finite duration 



32 



36 



40 



12 16 20 24 
time [ fm/c ] 

Fig. 2. The time dependent kinetic energy from a solution oi 
the Kadanoff-Baym equation (KB) is shown together with the 
results from the finite duration approximation ( |l3| ) and the 
Boltzmann equation. For Boltzmann transport the kinetic en- 
ergy is conserved in each collision and therefore globally. The 
conclusion is that the broadening of the 8 function of energy 
conservation in the finite duration approximation almost ac- 
counts for the time dependent built up of correlations from 
the exact KB equation. 



In figure g we compare the results with the exact solu- 
tion of Kadanoff-Baym equations pl|,p|. We see that the 
finite duration approximation reproduces the exact result 
quite nicely. The small deviation is due to higher order 
effects. 

In order to investigate a situation with higher degen- 
eracy we choose a model of two initially counter - flowing 
streams of nuclear matter with density and temperature 
ni = n o /60 Ti = 0.5 MeV and n 2 = n o /20 T 2 = 0.1 
MeV moving with relative velocity of lh/ fm correspond- 
ing to a collision energy of 21 MeV/n. The interaction is 
again a Gaufi-type of potential. In figure || is plotted the 
time evolution of the kinetic energy and the correlation. 

This build up of correlations is independent of the ini- 
tial distribution form. If for example we choose a (equi- 
librium) Fermi distribution as the initial distribution a 
build up of correlations will occur as well. This is due to 
the fact that the spatial correlations relate in momentum 
space to excitations, resulting in a distribution looking 
somewhat like a Fermi-distribution but with a tempera- 
ture higher than that of the initial uncorrelated Fermi- 
distribution. [^l| 

In order to illustrate the temperature dependence of r c 
as well as to demonstrate the quality of limiting analyti- 
cal formulae, we plot in figure ^ (thick lines) results from 
the solution of the Kadanoff and Baym equations for a 
fixed chemical potential of 37.1 MeV and for three differ- 
ent temperatures. The figure shows the increase of the ki- 
netic energy (equivalent to the decrease of correlation en- 
ergy) with time. The KB results are compared with those 
from approximation (|2^) . One sees that initially while cor- 
relations are built up the agreement is good. Especially at 
low temperatures the oscillations discussed above are ob- 
vious in the approximate results while the KB calculations 




10 20 30 

time [ fm/c ] 

Fig. 3. The time evolution of the correlation and kinetic energy 
for two initially counter -flowing streams of nuclear matter with 
density and temperature n\ — n o /60 Ti = 0.5 MeV and ri2 = 
n o /20 T2 = 0.1 MeV moving with relative velocity of lh/ fm 
which corresponds to a collision energy of 21 MeV/n. In this 
case the correlation time r c = 9 fm/c which is appreciably 
larger than in Fig We ascribe this to the smaller thermal 
velocity of the particles in the present case. The parameters 
are here in contrast to figure ji] chosen such that the system is 
degenerate nA 3 = 1.14 > 1. The equilibrated temperature is 
here (neglecting correlations) T = 11.6 MeV. 



only show a slight overshoot at the lowest temperatures. 
We believe that the discrepancy is due to the neglected 
damping (and perhaps due to the necessary approxima- 
tions used in the integrations as discussed above). The 
opposite approximate formula for high temperatures ( |l9| ) 
is plotted also for the T = 40 MeV case. The built up of 
correlations is too fast according this formula. 

From the numerical inspection of this section we sum- 
marize three facts: (i) the build up of correlations is mono- 
tonic and reaches the final value smoothly, (ii) the finite 
duration approximation ( |l3| ) is an excellent approxima- 
tion and (iii) also for initial Fermi distributions we see the 
same build up of correlations. The latter point shows that 
this build up of correlation is not due to an equilibration 
of an initial nonequilibrium distribution but mainly due 
to the decay of higher order correlation function which 
are condensed in the off-shell sin - factor in the collision 
integral (0). 



3.3 Correlation time 

From the numerical observation we conclude that the cor- 
relation energy (|l^) is increasing monotonously with time 
until it reaches its almost final value (|3^) . We assume that 
the main formation time of correlations r = t c — to is given 
by reaching this asymptotic limit. This time can be esti- 
mated by the condition 



-^corr (^~) 



1 ,COS 



AEt 



1 



4 AE 



AE 

F(oo)) 



-F(tc)) 
E 



(24) 
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Fig. 4. The formation of correlation plotted as an increase of 
the kinetic energy with time for temperatures 1, 10,40 MeV. 
The chemical potential is fixed to 37.1 MeV which corresponds 
to densities 0.16, 0.18, 0.35fm -3 . The thick lines show results 
from KB calculations while the thin lines are approximate val- 
ues via formula (jig]). The equilibrium correlation energy was 
chosen to be equal to the KB result. The oscillations are over- 
estimated by the approximate formula. For T = 40 MeV we 
plotted also the high temperature approximate value via (^) 
as thin solid line. The built up of correlations is too fast. 



We solve this equation for r approximately by replac- 
ing the cos-function by a linear approximation within the 
build up time interval t — to = (0, nh/AE) 



1 



AE(t-tp) 



AE 



2(t - tp) 



Sgn(AE). 



(25) 



This approximation is correct in three points, the initial 
time, the final formation time where < V(r) > has its 
maximum and a time point half of this maximal time at 
ttTl/2/AE. The latter point coincides too because 1— cosxt 
has there a turning point. Therefore this linear approxi- 
mation overestimates the function in the first half of the 
interval t — to = (0, ith/AE) and underestimates it in the 
second part. Furthermore we use equilibrium distribution 
functions. With the help of this approximation we can 
solve (§J) 



1 h 

2* (Sgn(AE)F(oc)) ' 



(26) 



have now 
27 



4T 



/ dxx 2 V 2 (x) J ds Sgn(x(s + x))c 

-oo 

OO OO 2 

fdxxV*{x) J ds^j 



(27) 



with V(x) = l/(x 2 + b 2 ) for Yukawa potential and V(x) = 
exp — (x/b) 2 for Gaufi potential. For the latter we obtain 
with the help of appendix ^ [(|62|)1 



TT 2 h 



(2 + 6 2 )(f + - arctan(\/2») 

fH { l- h l +0 ^)) 
8v t h 6 



(28) 



with the thermal velocity v 2 h — Tf[i and b = h/r)/s/2jjT. 

We see that in the low density or quasi-classical limit 
b — ► the formation time of correlations are determined 
entirely by the range of potential 7/ divided by the thermal 
velocity. This result is intuitively appealing. 

The opposite limit of low temperatures can also be 
evaluated like in appendix A. 2. This limit is independent 
of the used interaction because following (|4(]) the interac- 
tion part cancels out in (|2(]) . Then we end up with Fermi 
integrals which are similar to ( |68|) and the resulting for- 
mation time of correlations is via (E6|) 



TlowT 



3e/ + e c 



1+ 2 



(29) 



This time agrees with the time where the correlation en- 
ergy (E3I) has reached its first maximum 



2h 

e + e c 



(30) 



This correlation time limits the validity of quasiparti- 
cle picture which is established at times greater than this 
0. Incidentally, in the early 1950s the criterion H/hsT < 
t was supposed to limit the validity of the Landau Fermi- 
liquid theory for metals ]22] ]. Later it was shown by Lan- 
dau that this criterion is irrelevant and he proposed the 
correct criterion r > h/tp- Here we have explicitly cal- 
culated the formation time of correlations which is found 
to be equivalent to the memory time. It is to remark that 
this result describes just the memory- or collision duration 
time T mem = 1/E, see |T6| |, For nuclear matter at satura- 
tion density this time is about 4—5 fm/c and agrees with 
the numerical result of memory time ||23| . 



If we use < Sgn(AE)F >« £J corr < AE >, where < 
AE > is the mean transition energy of the collision, we 
obtain r w K ^ E> ■ This gives the intuitive picture of an 
uncertainty principle, i.e. a smallest time scale determined 
by h divided by the transition energy. 

The high temperatur e val ue of (|2^) can be calculated 
analogously to appendix A.l . Instead of (^) and 



3.4 Range of validity 

The validity of the low and high temperature expressions 
can be discussed with the help of two parameters, the 
value of the degeneracy nA 3 = cxp fi/T and the ratio e//T. 
There are 4 cases 
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Fig. 5. The 4 different areas of density- temperature range 
according to (|3l|). The high temperature expansion is applica- 
ble in case 2 and the low temperature expansion in case 4. The 
right figure represents the situation for spin-isospin degeneracy 
s — 4 and the left figure the situation for s — 1 for comparison. 



case 1 nA 3 > 1 T > et only for s > 1 

case 2 nX 3 < 1 T > tf 

case 3 nX 3 < 1 T < ef only for s = 1 

case 4 nX 3 > 1 T < e/ 

where s is the degeneracy. In figure (|[) we plot this case 
for nuclear matter with m = 938 MeV. 

It is clear that the low temperature expansion can only 
hold for case 4 and the high temperature expansion in case 
2. For nuclear matter with s > 1 we see that the case 1 is 
not covered at all by these expansions. 



3.5 Discussion 

We compare in Table @ the different expansions with each 
other as well as with exact results from KB-calculations. 
The general agreement of the numerical value of the cor- 
relation energy (132) with the KB result is striking. The 
correlation time (pfl) calculated numerically agrees as well 
reasonably. 

For the case 2 in row 1 (for a high temperature and 
nondegenerate system) we can reproduce the build up 
time quite well with the approximate formula ( p8| ) . This 
becomes worse if we approach the degenerated condition 
nA 3 — 1 in second line of the table. Both of the high tem- 
perature cases are too narrow to this condition to repro- 
duce the correlation energy correctly. The high temper- 
ature approximate formula ( |39| ) overestimates the exact 
value B C orr obtained from the solution of the time depen- 
dent kinetic equation. However in the case 4 for low tem- 
perature and degenerated systems we can reproduce the 
simulation value quite well. The build up time is slightly 
overestimated by the exact low temperature value (|30|). 
The case 1 is not covered by the approximate values which 
is demonstrated in the 5. row of the table. There the cor- 
relation energy is much overestimated which is quite clear 
from the discussion of Fig ||. The systematic slight in- 
crease of the correlation time with temperature (for low 
temperature limit) is explained by the low temperature 
expansion. 

The results related to figure |J are represented in rows 
6-8. It underlines the general good description of case 4 



and the failure of case 1. The chemical potential was kept 
constant initially in these three cases and the tempera- 
ture gradually increases, such that the chemical potential 
decreases. According to (|3fj| ) the correlation time should 
increase if using the final chemical potential instead of 
the initial one. This is clearly not the case in K/B so- 
lutions, where the correlation time stays almost constant. 
This shows clearly that the formation time cannot be com- 
pletely described by equilibrium formulae as done above. 
The correct formation time is described by a chemical po- 
tential somewhat between initial and final value. 

We would like to repeat here that the expression ( |2^ ) 
for the low temperature limit is universal for binary colli- 
sion approximations, independent of the interaction. The 
high temperature limit should also be correct because the 
Born approximation used here is believed to be a good 
approximation for fast particles. The intermediate region 
is quite open. We have calculated only in Born approxi- 
mation. Here especially higher order correlations beyond 
second Born should be employed. In summary: For low 
temperatures approximate analytical expressions can be 
given and for higher temperatures the analytical second 
Born approximation should be approached, while the in- 
termediate region is left for numerical investigations. 



4 Summary 

The gradient approximation of the kinetic equation in sec- 
ond order Born approximation is investigated. The inter- 
action energy is derived within different expansions. The 
equilibrium value of the correlation energy is obtained 
from the first order gradient expansion. This equilibrium 
value is calculated for Yukawa and Gaufi type of poten- 
tials and the results are analytically given in high and low 
temperature limit. For contact potentials we rederive the 
known result for the ground state correlation energy. 

A finite duration approximation of the non Markovian 
collision integral is proposed which follows from time de- 
pendent Fermi's Golden Rule and which is in good agree- 
ment with the numerical solution of complete collision 
integral. Furthermore it leads to the correct equilibrium 
value. Numerical comparisons are made with the solution 
of the complete Kadanoff and Baym equation in Born ap- 
proximation and with this finite duration approximation. 

The build up time of correlations is investigated and 
it is found that the low temperature value is universal 
for any approximation at the binary collision level. It is 
shown that the formation time of correlations is nearly 
determined by a ratio of h to the transfer energy which 
can be considered as an analogue to uncertainty principle. 

The high temperature limit shows roughly the time 
scale that particle needs to travel through the potential 
range. The validity of both the high and the low temper- 
ature estimates are confirmed by numerical comparisons 
with KB-results. 

The time scale we are describing are just the life time 
of fireball assumed in the early stage of nuclear collision. 
This means the extracted temperatures from final stage 
products are usually wrongly extrapolated to early stages 
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of nuclear reactions. This should effect the conclusions to- 
wards the caloric curve which is much discussed recently. 
Also the size and lifetime of hot reaction centers, which are 
extracted by interferometry methods should be critically 
revised as demonstrated in pij . 

The authors like to thank P. Lipavsky and V. Spicka 
for interesting discussions and A. Sedrakian for helpful 
comments. 



A Equilibrium correlation energy 

From the gradient expansion (]?]) we see that an equilib- 
rium value of the correlation energy =< VoF > is 
approached for large times. Here F is given by Eq. (||) 
with equilibrium distribution functions and Vq is given by 
( f43| ) such that we have 

EZAoo) - EZM = < ~Fo > 



sis 2 f dp 1 dp 2 dp[dp' 2 



niPi-p'il) 2 



h J (2tt^) 9 
x 5(pi +p 2 -p[ -P2)/o(p1)/o(p' 2 )/o(pi)/o(P2) 



Ei + E2 — E[ — E' 2 



(32) 



where in the second equality we have employed symmetry 
relations resulting in a factor 2 and have written all terms 
explicitly. 

We are now able to calculate the correlation energy 
analytically. In order to compare with known results we 
calculate first the ground state energy for contact poten- 
tials. Then we will be able to find analytical expressions 
for two extreme cases, i.e. the high temperature and low 
temperature limit. 

For performing the integrations in the collision inte- 
gral we define new coordinates by p = p[ — pi, II = 
T ^ 7 ( Pl + p'J and p' = p 2 p> 2 , II' = IS ^_ J ( P2 + 

p 2 ). For the 77, II' coordinates we use cylindrical coor- 
dinates II = Hp + p where U p is parallel to the z-axis 
fixed by p. Within these coordinates the momentum con- 
servation takes the simple form p' = p and the energy 
difference is AE = j^{n' p — H P ) with [i the reduced mass. 



For the contact potential we find from (32) the known 
ground state correlation energy |^5f 



E c 



3 4(21og2- 11) p f 



N 



2l7T 2 



<T> 2 



(33) 



with the scattering length a related to the coupling con- 
stant of the contact potential 4irh 2 a/m. As pointed out 
in p6| we have to subtract an infinite value, i.e. the term 
proportional to /i/ 2 in (|32]). This can be understood as 
renormalization of the contact potential and is formally 
hidden in 75^. (0) on the left hand side. For finite range 
potentials we have an intrinsic cut-off due to range of inter- 
action and such problems do not occur. For the standard 
procedure at low temperatures when the angular and en- 
ergy integrals are separated, we have to care about this 



cut-off once more because during this procedure we ne- 
glect the effect of potentials but we will have to integrate 
off-shell. 



A.l High temperature limit 

The high temperature limit is performed using the Maxwell- 
Boltzmann distribution ( ^8|) for the calculation of the cor- 
relation energy $2\\). 



A. 1.1 Yukawa-type potential 

Performing a series of trivial integrals we obtain the ex- 
pression for the correlation density for the Yukawa poten- 
tial 



El 



h[b] 



-*£(- 



4/i 



2b 



1 + mb 
xdx 



KT 

ds 



h[b] 



r 3/2 J (j.2 




b 2 f 



(34) 



where the abbreviation b 2 = (Tin) 2 /8/iT has been used. 
The parameter b proportional to h is responsible for quan- 
tum corrections. The last integral in ( |34| ) is performed in 
appendix ^ and the final result reads 



El 



4/x 



in,. 



) 2 ^p(l-V^be b2 eric(b) 



(35) 



This result is the correlation energy in second Born ap- 
proximation. The quantum corrections are condensed in 
the b dependent function. For identical plasma particles 
(g = e 2 /eo) the result reads 



E Y 



n 



e 

4e 



- fl - V^e b2 erfc(6) 
n V 



= ^l(l-^b + 2b 2 + o(h 3 )). 
4eo 



(36) 



The classical limit shows half of the exact low density 
limit of correlation or self energy, J27|, page 115, which we 
calculate in (B91) in appendix O 



^ e = ^^(l-e &2 erfc(6) 



2e 

— (1 - ^b +-b 2 
2e V 2 3 



+ o(h 3 )). (37) 



This difference follows from the static approximation 
of the screened interaction by the Debye potential, W(u) 
W(0) = Vq. The energy dependency of true W reflects 
that the screening of Coulomb potentials is a dynamical 
process. The Debye potential is a correct approximation 
for the scattering rates, since the transferred energy to is 
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small being limited by the on-shell condition Ae = 0. For 
the off-shell processes, the transfered energy can reach ar- 
bitrary high values for which the Debye potential would 
fail. Instead we would have to use the dynamical screened 
potential W(q, 2^ ) m the Levinson equation At 
low energies, the imaginary part of W needed in the scat- 
tering integral reads lmW = |W| 2 Imi7 w V^ImlT. The 
principle value integration in the energy function ( |32|) can 
be evaluated via Kramers-Kronig relations between the 
real and the imaginary part of the screened potential. 
The conjugated real part is Re(M / — Vc) = VfcReilW w 
VcRelJReW « VcRc77Vd, where 77 is the susceptibil- 
ity. Thus the dynamical screened potential would lead to 



, while the static approximation yields 
The factor of two in the exponent brings 



V C V ~ e- 

y2 _ e -2 K r T 

the factor of two into denominator in formula (1351). The 



corresponding integrals leading to (37) can be found in 
appendix |c[ This means that we do not consider the for- 
mation of effective Debye potential, but use it here as an 
effective one present from the beginning. 



the Yukawa- or GauB- potential reads 



E, 



( arctan ± 



h 

l+b'f 



Yuk. 



2V2n 



erf( 



El 



) Gaufl 



(42) 



where e 



/ 



e f +e c 



and bj 



1 u l = ip' ^ es * cn °i ce f° r cut_ 

off we found e c ~ ej. This result shows that even for zero 
temperature we obtain a correlation energy due to ground 
state correlations by Pauli blocking, which we have explic- 
itly demonstrated for contact potentials ([33j). 

It should be noted however that the temperature T 
refers to the final equilibrated distribution. It is, in gen- 
eral, not a Fermi-distribution but is modified because of 
the correlations. This is for example illustrated by the re- 
sults shown in figures 1 and 3 in pi"j. 



A. 1.2 GauB-type potential 



A. 3 Numerical values 



For the GauB potential the same steps as before lead in- 
stead of (|U) to the expression 

OO OG 

23/2 f X e-^ dx ( -*L e -' a , (38) 



h[b] = |p- / xe- 2 ^dx J 



s + x 



where b 2 = n 2 /{2^Trj 2 ). With the help of (||) in appendix 
|^ we obtain 

F G ( 4 M s 2 ^ /2 Von i n a n b 1 

^ma + mj 2V2T 2 + 6 2 ' 1 ' 



A. 2 Low temperature limit 

In the opposite limit of small temperatures compared to 
the Fermi energy we employ the well known methods of 
Fermi liquid theory p8| , p9| . For simplicity we now use only 
one type of particles. The angular and energy integrals in 
(p2\) can then be separated with the result 



1 V 2 

E corT = -~m Pf( W )If, 



(40) 



where 



V 2 , 1 [d fhdfMfo 2 

' a — v ypi-Pi) (4i) 



(27T^) £ 



with standard notation of angular integrals fl28f . The Fermi 
integral If is done in appendix [5] and the final result for 



In order to check these low and high temperature ap- 
proximations, we plot in figure I the results from @ 
and ([39]) with the numerical result of ([32]) for the GauB- 
potential. We have chosen a normal nuclear density, n = 
0.18 fm~ 3 . We see that the low temperature expansion, 
(ff2|), gives reasonable agreement with ( |32|) at low tem- 
peratures only. The high temperature limit is approached 
actually only at very high temperatures beyond 100 MeV 
at this density. We would like to emphasize that the nu- 
merical computation of of (|3|) involves a complicated four 
- dimensional principal value integration. While we con- 
sider the high temperature limit to be accurate, the low 
temperature limit is approximate and dependent on the 
energy integration. The cut-off value of the energy inte- 
gration (see discussion after (|32])) was chosen to be e/ 
throughout the calculation. In the low temperature limit 
we used instead of e/ the Sommerfeld expansion fi(T) = 

e /( 1 -S($)V°c r8 )- 



B Weakening of initial correlation and the 
limit of complete collisions 

In order to derive kinetic equations for the long time evo- 
lution one uses the limit to —> ~~ 00 where to is the initial 
time where correlations beyond the one-body are assumed 
to be absent. Two- particle correlations vanish then in the 
remote past. ^,From (U) it is obvious that we could con- 
sider t — > 00 alternatively. This is sometimes referred as 
complete collisions , i.e. uncorrelated asymptotic states, 
which can include one-body correlations, of course . 

We like to show now that special care is needed when 
this limit is carried out, especially, if one likes to consider 
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Fig. 6. The correlation energy (|32|), solid line, is compared 
with the low (|l2|), dotted line, and the high temperature (^), 
dot-dashed line, approximations. Results are shown as a func- 
tion of temperature at a nuclear density of n = 0.18fm -3 . See 
also discussion in text. 



non-markovian corrections as we did within the gradient 
approximation (Q) . Especially we will find that the order- 
ing of the gradient expansion and this long time limit are 
not easily interchangeable. 

To illustrate this fact let us first carry out the to - * 
— oo limit on the zero and first order terms V n (t) in Eq. 
(0). This implies that we consider the limit of complete 
collisions. The result is 



AE(t-t ) 



V Q (t) = 

Vi(t) = 



AAE 
h d sin 



AE(t-to) 



1 P 

4 AE' 
h 



4 d AE 



AE 



S'(AE), 



V n (t) oc h n , 



(43) 



where P is the principal value as given in the appendix. 

Therefore we obtain in the long time limit (i.e. for 
completed collisions) 



d 

— lim i? CO rr 



(44) 



On the other hand we can first perform the time derivative 
and then the long time limit. To this end we employ the 
gradient expansion of the collision integral (^|) and obtain 
with the help of (0) 



d_ 



£■„„, = < I F(t) > + < I x F'{i) > +o{F"{t)) 

o AE(t-t ) 



18_ 

45t 



< 



1 



to) sin 



AE 

. AE(t-t ) 



> 



F(t) 
F'(t) > 



9 1 P s 1 d P , 

--0-t < lAE F ^ > - < r07AE F V > > 



(45) 



where the relations T cos xT — ► — e-^nS e (x) andTsinxT — » 
— e ■§-— has been used. We see that we obtain with the help 



of x 



de x 
P' 



lim 



dt 



-JLE 

-de x 



Ecorr = \< AE^F>(t) > +0(F"(t)), 



(46) 



where P' is the derivative of the principal value as given 
in the appendix. Comparing with (^) we recognize that 

d d 

t — »—oo Ot (Jtto—>-°a 

-<4^*«»- 

(47) 

Here we can deduce the rule from ( pT| ) for large times 

£grad-c*act W >= tfg-d-coU (j. ) + < ^ £ 9__£_ F(f) ^ 

that the exact correlation energy in gradient approxima- 
tion for large times i5|™ r d_cxact (t) is given by the corre- 
lation energy from gradient expansion of collision inte- 
gral ES^~ coU (t) via an infinitesimal correction term. For 
any function F with existing derivatives the distribution 
e ~cTe~3E~ vanishes- Therefore we are allowed to identify the 
correlation energy from the kinetic equation in gradient 
approximation with the corresponding correlation energy 
in second Born approximation. Total energy is conserved 
in first order gradient expansion for nondegenerate sys- 
tems [pl| . For the degenerate systems we can prove com- 
plete energy conservation also for degenerate and space 
nonlocal systems . 

However we like to point out that the above result is 
not valid if the transition probability contains poles, e.g. 
bound states or pairing ones. Then we get a nonvanishing 
last term in ( |48| ) expressing long living correlations such 
as bound states or pairing. The same discussion is applica- 
ble for formal divergent expressions as it appears e.g. for 
contact potentials. With the help of ( f48| ) we can renormal- 
ize the expression to obtain the ground state energy as we 
will demonstrate later. 



B.l Long time limit of the collision integral 

The message of the last paragraph does not mean that 
the long time limit to — > — oo of the collision integral is 
not uniquely defined. To show this we consider the two or- 
derings of operations of the time derivatives and the limit 
of remote past respectively. For this purpose we write the 
gradient expansion up to first order of the kinetic equation 
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in two different ways 

F(t) 



2b 



~ , ■ (t-t )AE 

d_ _ 1 sin I 
dt h ~ h AE 



cdx 



R J (x 2 + b 2 ) 2 
o 



lmw(—x + ie), (53) 



d cos 



AE(t-t ) 

h 



dAE AE 

, / • (t-t )AE 
1 / sm o> 



-F'(t) + o(F"(t)) 



AE 



(t-to) AE(t-t ) 
— - cos 



d 

+ dt 



o AE(t-t a ) -, 

a cos — it — 11 - 1 



dAE 



AE 



F(t) 



F(t) 
(49) 



The first equality is written in terms of the derivatives F 71 
Now we take the limit to — > — oo on this and obtain 



8t 



1 



he 2 + (AE) 
J*«US)F<*>-| 



at 



8 



AE 



dAE (AE) 2 + e- 



F(t) 



P' 



(50) 



Equation (pO) is just the generalized kinetic equation de- 
rived in |16|,|17[ which leads to quantum Beth-Uhlenbeck 
virial corrections in nonequilibrium. It consists of the usual 
Boltzmann collision integral and a correction represented 
by the second term on the right hand side of Eq. (50). 

We next invert the order of operations, i.e. the time 
derivation and the long time limit, and therefore use the 
second equality of (^) where the time derivative is in front 
of the collision integral. Applying the limit to — * — oo one 
obtains 



±f =-(- 

dt n h { e 2 



0_ 

Iff 



d 



f- (AE) 2 de e 2 
d AE 



(AEf 



)F(t) 



dAE (AEf 



F(t) 



(51) 



Whereas the second term is just the same principal value 
expression as in (|50|), the first term is at first glance dif- 
ferent, but one finds 



8 



2e 2 



e 2 + (AE) 2 de e 2 + (AE) 2 e 2 + (AE) 2 e 2 + (AE) 2 

= ttS(AE) (52) 

so that the two results (^0|) and ( pl[ ) are in fact equivalent. 

This means that despite the care which is necessary for 
obtaining the mean value of the correlation energy (fl8|), 
the first order gradient expansion for the generalized ki- 
netic equation has a unique long time limit independent 
of the ordering of the two operations involved. 



C Error function integrals 



The following integral is needed in (35) for the calculation 
of the correlation energy 



where the complex error function has been introduced [p0| 



w(z + te) = — 



7T J z — t + ie 



(54) 



Performing a partial integration and using w'(z) — —2zw(z) 
2?/y / 7r we get 

oo 

2b f xdx 
h[b) = l + — / - Imw(-x + ie). (55) 

Vtt J (x z + b z ) 
o 

Reintroducing (|54|) the last integral is easily carried out 
by interchanging the two integrals and we obtain finally 

h[b] = l-6\A r c fc2 crfc(6). (56) 
In order to verify ( |37| ) we have instead of (|5J 



7T 3 / 2 J X(X 2 +b 2 ) J S + 











OO 


dx 


r ds 


X 


J s + X 




OO 


ds 




+ b 2 





(57) 



where a simple decomposition has been performed. For 
the first part we apply the Poincare- Bertrand theorem 
[Bill which states 



dx f dy 



x — u J y — x 
dy / dx 



■f(x,y) 



f(r ' ,n ir 2 f(u,u) (58) 



(x - u)(y - x) 



and obtain 



hi = ^y(l-e b erfc(6)). 



(59) 



Next we perform the integral required in (|38|) where 
we put c = 2/6 2 



ds 



s — X 



h[b] 



xdx 



2b 

^T 2 J (x 2 + b 2 ) 2 
o 



ds 
( 

s + x 



= _ _L / ^ e -* s2 (/( S )-/(- S )) 

J s 

o 

oo 

f(s) = f dxxe- {1+c)(x+ ^ )2 . (60) 
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Here we have derived the last form by dividing the integral 
s into two parts corresponding to the pole at s = ±x 
and performing variable substitution. The integral f(s) is 
elementary and yield 



-erfc(- 



1 



(1 + c) 3 / 2 2 V % /TT7 ; 2(1 + c) 



e <!+<=) . 



(61) 



Substituting this into we obtain the result 

on 

2 3/2 



h[b] = 



6 3 (l + c) 3 /20f 
V2 1 



e !+<= s ds 



b 3 y/c 1 + C 

1 



6 2 + 2' 

In (Kfl) we need the integral 



(62) 



7 4 [6] 



46 



x dx 



^3/2 y ^2+^2)2 / ^ " Sin( S + x)xi 



-t 2 a: 2 /4 



26 


oo 

/ 


7T 






OO 


1 


d 


7T 


db 


Im 


(1 



(x 2 + b 2 ) 



-(t 2 /4+it)x 2 



(l + 2z 2 )c z (l-erf(z)) 



2z 



(63) 



where we used z — byjt 2 /4 — ii. 

Another error function integral is used in ( |28| ) 

OO 



(64) 



To this end we use a tabulated integral Q (6.285.1) 



la! [a] = / dxe" * erf (a:) 
o 

? — arctan-^/a 



(65) 



The desired integral is then 



h[a] = -Vn-z-hi[a] 
a a 



1 /tt 



2a 3 / 2 V 2 



arctan(va) 



(1 + a) 



(66) 



D Calculation of collision integral 

In order to calculate the collision integral for Fermi func- 
tions in ( |40| ) special care is needed. The problem differs 
from the usual calculations of collision integrals in that we 
do not have energy conservation but an off-shell principal 
value. Therefore the limit A = fi/T — > oo cannot be per- 
formed straightforwardly. The reason is that the energies 
are restricted to the neighborhood of the Fermi level, but 
we have to allow off-shellness due to the denominator of 
the energy expression. We have used the approximation 
to separate the collision integral into angular and energy 
integrations. This leads to formally divergent results be- 
cause the matrix element of interaction is not included in 
the energy integrals and the latter become infinite. In or- 
der to cure this defect we restrict the energies available 
from above. Then the expressions become analytical and 
we fix the upper limit of energy integration by match- 
ing to the numerical result found from (|32|). The best fit 
was found to be two times Fermi energy, i.e. A c = A. We 
shortly sketch the derivation with this upper limit. 

Introducing dimensional variables x — (e — /i)/T and 
A = fi/T we have 



If = T 3 



dx\dx2dxdx' 2 



(e*i + l)(e*= + l)(c- x 'i + l){c- x 2 + 1) 



(xi + x 2 - x[ - x' 2 ) ' 



(67) 



Now we switch to difference and center of mass coordinates 
via p = (x\ + x 2 )/2 and r = x\ — x 2 and analogously for 
p',r'. The integrals over r and r' are then trivially carried 
out because their integration range is (— (A+A c ), A+A c ) 
(—00,00). In new coordinates x = p 
the result reads 



p' and z = p + p' 



A+A c 2X r 

-e~ x I dz- 

X 

-(A+A e ) 



h = -\ T " ! -*- x I dz f' z \ 

2 I x I cosh x — cosh z 



-2A 



2(A+A C ) 
-2(A+A C ) 

47r 2 r 3 (A + +A e ) + 1. A + A e a 

3 3 7T 



(68) 



where we have employed the integral ([73j). We see the ex- 
plicit dependence of this off-shell integral on T 3 A = T 2 ef. 

We would like to remark that the limitation of up- 
per available energies is equivalent to cut- off procedures, 
which has been employed within memory collision inte- 
grals. This cut- off comes here from the approximate cal- 
culation of the collision integral and is fixed to reproduce 
the numerical results. 
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Fig. 7. The integration contour for the integral ( |73| ) 

E A useful integral 

To evaluate the useful integral 

oo 

J„ = P / dx 



cosh x — cosh ( 



(69) 



we use the integration contour as depicted in figure 0. 

We have the following identity in the limit of large R 
and small radius around the poles 



, , x n - (x + 2iri) n , , 

P / dx u J =- Cl -c 2 ~c 3 ~c A (70) 

cosn x — cosn a 



where the Cj are the infinitesimal half circles surrounding 
the poles at ±a. The latter ones are trivially calculated 
with the final result 

oo 

P / dx Xn - {x + 2lTl)n =J^(a n - (-a)" 
J cosh x — cosh a sinh a 

— oo 

+ (2ni + a) n - {2m - a) n ) . 

(71) 

From this identity one now generates the required inte- 
grals J„ by calculating the identity ([7l] ) for n + 1. The 
first two integrals read 



Jo 



J-2 



2a 

sinh a 
2a 

3 sinh i 



-(a 2 -27r 2 ). 



(72) 



As is obvious from symmetries only integrals J n with even 
numbers of n are nonzero. The required integral in ( p8[ ) 
reads now 



P / dx- 



cosh x — cosh a 



= J 2 - a 2 J 



4a 



3 sinh a 



(tt 2 + a 2 ). (73) 
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Table 1. Comparison between estimates discussed in text and KB- results. The first two rows are calculated with the situation 
of figure |l]and Fig|. The last six rows refer to results with Fermi distributions as initial conditions. In the first three of these the 
density was kept fixed and the temperature was changed. In the last three the chemical potential was fixed and the temperature 
was changed. Temperatures are given before and after the formation of correlations. The latter temperatures are higher due to 
the (negative) correlations. For the correlation energy _B C orr and the correlation time r c we compare 3 different expressions with 
the K/B solution. The numerical values (^) and(pi|) has been performed by four-dimensional principal value integration. The 
approximate formulae for the low temperature limit (W2J) and (pOj) are given as well as the high temperature limit (pflj) and (Eq). 
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